The Polytropic Equation of State of Primordial Gas Clouds 
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ABSTRACT 

The polytropic equation of state (EOS) of primordial gas clouds with modest emichment is 
computed, motivated by the recent observations of very Fe-deficient stars, [Fe/H]^ 10^^'^ — 10~^, 
such as HE 0107-5240 and CS 29498-043. These stars are overabundant, relative to Fe, in C and 
O. We assume that the observed abundances of species like C, O, Si and Fe are representative of 
the gas from which the currently observed metal-deficient stars formed. Under this assumption, 
we find that this primordial metal abundance pattern has profound consequences for the thermal 
balance and chemical composition of the gas, and hence for the EOS of the parental cloud. The 
polytopic EOS is soft for low, [0/H]< lO^'^, oxygen abundances, but stiffens to a polytropic 
index 7 large than unity for [0/H]> lO^'^ due to the large opacity in the CO and H2O cooling 
lines. It is further found that a regulating role is played by the presence and temperature of the 
dust, even when the overall carbon abundance is only [C/H]~ 10~^. When the dust is warmer 
than the gas, a region with 7 ~ 1.2 results around a density of ^ 10'* cm"'^. When the dust is 
colder than the gas, a region with 7 ^ 0.8 is found for a density of ~ 10^ cm"'^. Implications for 
the primordial initial mass function (IMF) as well as the IMF in starburst galaxies, where the 
metallicity is super-solar, are explored and related to processes that influence the temperature of 
the ambient dust. 
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A fundamental issue in the study of star for- 
mation is understanding the physical structure of 
molecular gas clouds from which stars are formed. 
Despite a wealth of observational data for the 
Milky Way (c.f. Fuller & Myers 1992; Goodman 
et al. 1998), the nature of the equation of state 
(EOS) remains a major theoretical problem in un- 
derstanding the stabillity and collapse of molecu- 
lar clouds. The stiffness of the EOS can be largely 
responsible for the resulting density probability 
function of interstellar gas in the turbulent ISM. 
In particular, the value of the polytropic index 7 
strongly influences the mass distribution of den- 
sity condensations as well as the amount of clump 
fragmentation (Vazquez-Semadeni et al. 1996; Li 
et al. 2003). In Spaans & Silk (2000) the proper- 



ties of a polytropic EOS were investigated and it 
was found that the stiffness of the EOS depends 
strongly on the ambient metallicity. 

Recent observations of metal-deficient stars 
such as HE 0107-5240 (Christheb et al. 2004) and 
CS 29498-043 (Aoki et al. 2004) indicate that 
formation of (low-mass) stars is possible in very 
low metallity gas clouds. The formation of stars in 
low-metallicity environments (Abel, Bryan & Nor- 
man, 2002) is important for the reionization of the 
universe and the early enrichment of the interstel- 
lar and intergalactic medium (Hirashita & Ferrara 
2002). Even though more (but not many more) 
stars are known with very low Fe abundances (e.g., 
CS 22949-037; Aoki et al. 2004), we concentrate 
on the two stars mentioned above since they ap- 
pear to reflect the wide range in C, N, O and a 
element abundance patterns for Fe-deficient stars. 
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It is not generally known what influence the 
environment has on the shape of the IMF. 
HST/WFPC2 observations of the LMC by Goulier- 
mis, Brandcr & Henning (2004) indicate that local 
conditions seem to favor the formation of higher 
mass stars (top-heavy IMF) in associations, and 
not in the background field. However, Meyer et al. 
(2004, and references therein) find that the IMF 
is not a strong function of environment or initial 
conditions. In any case, the numerical simulations 
of Li et al. (2003, and references therein) indicate 
that the value of the polytropic index strongly 
influences the spectrum of density condensations 
that is formed under the influence of gravity and 
turbulent driving. In particular, fragmentation 
is enhanced or inhibited when 7 is smaller (0.7) 
or larger (1.1) than unity. In this context, we 
present computations of the polytropic EOS of 
primordial (low-metallicity) gas clouds at high 
redshifts. We also comment on interstellar re- 
gions that have super-solar metallicities and are 
exposed to a warm (dusty) infrared background 
(e.g., starbursts). 

2. Model Description 

We assume that the observed abundances of 
species like C, O, Si and Fe are representative of 
the gas from which the currently observed metal- 
deficient stars formed. We further employ a poly- 
tropic EOS, P GC p^, where 7 is the polytropic 
index and p the mass density. We adopt a per- 
fect gas equation of state, P oc pTg for the gas 
temperature Tg, to write 7 as 

dlog7^_ 
cflogp 

This last step is justified (Scalo & Biswas 2002; 
Vazquez-Semadeni, Passot & Pouquet. 1996) as 
long as the heating and cooling terms in the fluid 
energy equation can adjust to balance each other 
on a time-scale shorter than the time-scale of the 
gas dynamics (i.e., thermal equilibrium). It should 
be noted that because 7 depends on the (loga- 
rithmic;) derivative of the temperature with re- 
spect to density, it implicitly depends on radiative 
transfer effects and changes in chemical composi- 
tion through derivatives of the heating and cooling 
functions (Spaans & Silk 2000). The model de- 
scribed in Spaans & Silk (2000) is used in this work 
and the interested reader is referred to that paper 



for a detailed description of the various heating 
and cooling terms that influence the polytropic 
index. The (updated) two-dimensional numeri- 
cal radiative transfer code of Spaans (1996) with 
the detailed chemical treatment of Spaans & van 
Dishoeck (1997), and its extensions to lower metal- 
licities by Norman & Spaans (1997) and Spaans & 
Norman (1997) is used in this work. With the 
work of Scalo (1998) and Le Bourlot et al. (1999) 
it is possible to extend the Spaans & Silk (2000) 
model to metallicities much smaller than 1% of 
solar, where H2 and HD completely dominate the 
cooling of the gas. The main features of the model 
are summarized below. 

We consider self-gravitating spherical clouds 
and therefore we adopt a singular isothermal 
sphere (Neufeld, Lepp & Melnick 1995), for which 
the total hydrogen number density nn scales with 
column density A^h per unit of velocity as 

Nu = 7.2 X 10^^n^^cm-'^/{kms-^). 

We include number densities, rin = n{H) + 
2n(H2), upto 2 X 10^ cm^'' and metallicities as 
low as 3 X 10"** and in solar units. Furthermore, 
dust temperatures between 10 and 100 K can be 
treated and the latest chemical reaction rates have 
been used in the computations (Le Teuff, Millar 
& Markwick, 2000). Care has been taken to check 
the convergence of the chemical equilibrium cal- 
culations for 0/C ratios that significantly differ 
from solar. In general, the thermal and chemical 
balance was solved iteratively with a convergence 
criterion of 0.1% for the chemical abundances and 
level populations. 

The computed grid is accurate enough to do 
linear interpolations between the log of any two 
adjacent values. Extrapolation of the values of 7 
to densities lower than 10^ cm~^ (say to 1 cm~^) 
is reliable because the slopes are modest. Note 
that softening of the EOS due to gas cooling on 
cold dust can be important and flips over to gas 
heating when the dust is warm. 

We have used a cosmic ray heating rate for the 
low-metallicity models that is given by a cosmic 
ray ionization rate of with ( = 3 x 10~^^ s~^. 
This means that cosmic ray heating is provided by 
a background star formation rate of the order of 
IMqIjv. However, the zero metallicity polytropic 
index is only weakly dependent on temperature 
between 100 and 2000 K so that the actual value 
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of the cosmic ray heating rate, as long as it is pro- 
portional to density, does not matter much for the 
low-metallicity case (Scalo & Biswas 2002). For 
starburst environments, we assume a cosmic ray 
ionization rate of 30(. 

Solar relative abundances are assumed (As- 
plund et al., 2004; Jenkins 2004). The latter con- 
dition is replaced by the observed abundance pat- 
tern of metal-deficient stars for some of the pri- 
mordial gas cloud models. We further adopt a 
MRN grain size distribution (Mathis, Rumpl & 
Nordsieck 1977) and assume that the dust abun- 
dance scales with the carbon abundance. 

The velocity dispersion of microscopic turbu- 
lence is AVtur — 0.5 km/s in the quiescent primor- 
dial gas, to which the purely thermal contribution 
is added (in the square) while the iterations on 
the thermal and chemical balance are performed. 
However, for the starburst model AVtur = 3 km/s 
is adopted to take the larger input of kinetic en- 
ergy (e.g., through supernovae) into account. 

No freeze-out of molecules is assumed. The 
latter assumption is justified because our dust is 
generally warmer than 20 K. The total extinc- 
tion through the starburst region containing the 
model cloud is 100 mag (Spaans & Silk 2000). 
The results of Le Bourlot et al. (1999) are used 
for the collisional excitation of H2 by H, He and 
H2, where non-local thermodynamic equilibrium 
level populations and quantum-mechanical cross 
sections are computed. For HD cooling, the re- 
sults of Flower et al. (1999) arc used. Both H2 and 
HD line emission are optically thin. For complete- 
ness we have also included H2 formation heating 
with equal contributions from the binding energy 
(4.93 eV) to translational, rotational and (statis- 
tically distributed) vibrational degrees of freedom 
(Sternberg & Dalgarno 1989; Meijerink & Spaans 
2004). 

The computed values of 7 depend sensitively on 
a number of effects: 1) line trapping, 2) metallic- 
ity, 3) H2O heating and cooling, 4) microturbu- 
lence versus a large velocity gradient (LVG) veloc- 
ity field. In particular, 4) is important in the con- 
text of collapsing molecular cloud cores. Note that 
a LVG approximation suppresses line trapping 
somewhat, and yields a relatively softer equation 
of state (see Spaans & Silk 2000). Also, H2O heat- 
ing, where absorption of infrared photons and col- 
lisional de-excitation causes heating (Takahashi, 



Hollenbach & Silk 1983), plays a smaller role for 
a LVG model unless Ay > 20 mag. 

We assume a microturbulent velocity field for 
all simulations. For optically thin (e.g., low- 
metallicity) line emission the velocity field plays 
no role. We ignore time-dependent effects (Bromm 
& Loeb 2003). UV illumination from the outside, 
not included here, will create a warm FDR region 
where the cooling is dominated by [CII] and [01] 
fine-structure emission and the polytropic index 
varies between 0.5-1.0 for C+ and between 0.6-1.2 
for 0-dominated emission. 

3. Model Results 

In order to assess the impact of elemental abun- 
dances that strongly deviate from solar, a grid was 
run with the parameters listed in Table 1. The 
figure number and color coding of each model are 
listed. The primordial model 1 is used as a base- 
line to assess the impact of C, O and dust on the 
polytropic EOS. Redshifts of 10 and 20 are con- 
sidered for the high redshift universe and the tem- 
perature of the dust, if any is present, is taken 
equal to that of the CMB. Overall, the grid cov- 
ers the two Fe-deficient stars, through models A, 
Bl and B2, and a range in carbon that limits, 
through model 13, to the zero-metallicity primor- 
dial cooling curve. Furthermore, various oxygen 
models are explored to show the large impact that 
O has on the EOS, i.e., models 10-12 for no O 
and models 2-9 for solar and super-solar O/C ra- 
tios. The latter are relevant because observations 
of CS 29498-043 (Aoki et al. 2004) indicate that 
the oxygen over carbon ratio can also be enhanced 
relative to the solar ratio. Finally, model 14 is a 
starburst environment where O/C is solar, but the 
overall metallicity is super-solar by a factor of two. 

It is straightforward to identify a number of 
trends in the models. We first note that in all 
models with a non-zero oxygen abundance, part of 
the gas cooling is provided by the [01] 63 um line. 
This line can go into absorption against the warm 
CMB below a gas kinetic temperature of ~ 80 K. 
Collisional de-excitation of excited [01] then leads 
to heating and a change in the polytropic index 
with the CMB temperature, e.g., in models 2-5, 
even when no dust is present. The same effect 
occurs for absorption of the CMB by water, i.e., 
H2O collisional de-excitation heating. From fig- 
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Table 1: All models use AVtur = 0.5 km/s, except for model 14, which uses AVtur = 3 km/s. The observed 

O and C abundance for the halo stars HE 0107-5240 and CS 29498-043 arc adopted in models A, Bl and B2. 
For CS 29498-043 the two adopted C abundances are indicative of the observational uncertainties in carbon. 
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ures 1-4 one sees further that: i) Carbon alone, 
even at 3% of solar, is not able to significantly 
stiffen (7 > 1) the EOS (figure 3). ii) Oxygen 
needs to be present at no less than the 1% of solar 
level to boost the opacity in the CO and H2O lines 
at densities around 3 x 10"^ — 3 x 10** cm~'^, see fig- 
ure 1, and yield values of 7 larger than unity. Note 
in this that the total abundance of oxygen (rather 
than carbon) sets a limit to the maximum amount 
of gas phase CO when [0/H]< 1/2[C/H]. in) The 
model parental cloud of HE 0107-5240 is charac- 
terized by a soft EOS, while that of CS 29498-043 
has a stiff EOS, see figure 4. iv) An overabundance 
(relative to solar) of oxygen leads to a significantly 
stiffcr EOS (figure 2), even when the overall mctal- 
licity is only of the order of 3% of solar. That is, 
if the O/C ratio is enhanced by a factor of three 
compared to its solar value, then the EOS stiffens 
to values of 7 in the 1.2 range at around tzh ~ 10^ 
cm~^. v) This effect is further enhanced (also fig- 
ure 2) when the temperature of the dust is high, 
due to the CMB at z = 10 - 20, and allows H2O 
(collisional de-excitation) heating to dominate the 
thermal budget and cause the gas heating to attain 
a rig density dependence, vi) At densities around 
10^ cm~^, the polytropic index dips below unity 
if the dust temperature, T^, is lower than the gas 



temperature. This results in 7 ~ 0.8 for models 
with dust at 2; = 0. If the dust grains are warm, 
then dust-grain heating causes the EOS to stiffen 
further and 7 Ri 1.1 beyond the H2O opacity peak. 
Recall in this case that gas-dust heating/cooling 
scales as (Tg - Td)T^-^n\. vii) The starburst (fig- 
ure 4) has the stiffest EOS because of its high abso- 
lute metallicity and warm (T^ = 75 K) dust, viii) 
The temperature of the dust, if present, plays an 
important regulating role, see figure 2. 

4. Discussion 

In commenting on our results, we wish to dis- 
tinguish slightly enriched, Z/Zq ~ 10~^, primor- 
dial gas that leads to Fe-deficient stars (models 
2-12, A, Bl and B2) from almost zero metallic- 
ity, Z/Zq ~ 10~4, primordial gas associated with 
pop III star- forming regions (models 1 and 13). 
As figures 1-4 show, these two metallicity regimes 
exhibit quite different behaviors as far as the poly- 
tropic index is concerned. 

We point out that all the models presented in 
this work are computed for a microturbulent ve- 
locity field with AVtur = 0.5 (or 3 for the starburst 
model 14) km/s and a thermal velocity contribu- 
tion that is determined by the thermal and chemi- 
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Fig. 1. — Figure 1 shows models 1-5; note that 7 
stiffens to just above unity around 3000 cm~^ for 
an oxygen abundance of O.OIZ0. 




Fig. 2. — Figure 2 shows models 1 and 6-9; note 
that 7 is well above unity for densities above 3000 
cm~'^, driven by the oxygen abundance. 




Fig. 3. — Figure 3 shows models 1 and 10-13; note 
that 7 remains below unity for densities upto ~ 
10^ cm-3. 



cal balance of the medium. However, if large veloc- 
ity gradients are present then optical depth effects, 
particularly line trapping in H2O and CO transi- 
tions, are more modest and 7 decreases (Spaans & 
Silk 2000). 

Also, adiabatic (compressional) heating has a 
density dependence of Fa oc n^^, in between cos- 
mic ray ionization, nj^*', and gas-dust or H2O, n^'^, 
collisional heating. Compressional work is not rel- 
evant to model clouds that can be described by 
singular isothermal spheres; however, in a collaps- 
ing medium it must play a role. Our heating terms 
bracket adiabatic heating when dust is present 
and/or Z > 10'^ Zq. For Z/Zq < 3x 10~*, we ex- 
pect 7 to attain a slightly larger asymptotic value 
of 1.0 - 1.3 for riH > 10^ cm'^ (Scab & Biswas 
2002; Schneider et al., 2002). At smaller densities 
the higher resulting temperatures, Tg ~ 200 — 800 
K (Abel, Bryan & Norman, 2000), prevent 7 from 
increasing significantly under the Fa term because 
excited atomic (e.g., O, Si) and molecular (e.g., 
CO, H2O) states with progressively higher critical 
densities become accessible, already at a metallic- 
ity of 1Q~'^Zq. These higher critical densities en- 
sure that the sub-thermal density dependence 
is maintained for larger gas densities. 

4.1. Fe-Deficient Stars 

In Spaans & Silk (2000), it is argued that a 
defining criterion for low- or high-mass star for- 
mation to prevail, i.e., whether fragmentation can 
proceed to small enough mass scales to cause a 
flattening in the IMF below about 0.3 Mq (Scalo 
1998), is whether the EOS has a polytropic in- 
dex 7c that is smaller or larger than about unity, 
respectively. Above this canonical value of 7c = 
2(1 — 1/1^), with v the number of dimensions in 
which cloud contraction occurs, sheets [v = 2) are 
able to stop fragmentation and start gravitational 
contraction because the Jeans mass does not de- 
crease any further. In the present discussion, the 
thermal Jeans mass, Mj oc T^-^ /p'^-^, as well as its 
generalisations to turbulent and magnetic media, 
plays a secondary role. 

If the elemental abundance pattern that is ob- 
served for HE 0107-5240 is typical for the parental 
cloud in which this star formed, then low-mass 
star formation is the dominating mode. This is 
consistent with the modest, 0.7AIq, mass of HE 
0107-5240. Furthermore, for a solar 0/C ratio, 
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a carbon abundance of less than 3% of solar still 
results in an EOS where 7 does not become sig- 
nificantly larger than unity even for warm dust 
grains. This would imply that the IMF would not 
be top-heavy at high redshift once metal enrich- 
ment sets in. Such an IMF would not be desirable 
for scenarios where the bulk of the reionization is 
caused by stars. We also checked that the abun- 
dance pattern observed in HE 0107-5240 provides 
sufficient levels of carbon to ensure that the cool- 
ing time is less than the free-fall time (Bromm & 
Loeb 2003). 

However, if the abundance pattern of CS 29498- 
043 is typical, then then we would expect that 
oxygen is relatively overabundant and that the the 
EOS is stiff, particularly when the ambient dust is 
warm. This would cause a shallow IMF at high 
redshift and a continuation of the formation of 
high-mass stars. Such an IMF would yield a larger 
number of ionizing photons per unit of mass. 

In this context, it should be mentioned that the- 
oretical models for supernova yields of very mas- 
sive, 20 — I3OM0, primordial stars are capable of 
reproducing the C and O abundances, but at the 
expense of the r- and s-process elements (Umeda 
& Nomoto 2003). An alternative would be that 
massive stars with M < 8Mq, i.e., stars that do 
not go supernova, evolve rapidly enough to pro- 
vide He-burning elements without the usual su- 
pernova dust production (Pel, private communi- 
cation). The latter scenario is interesting in that 
one would have a stiff EOS around ^ 10"* cm~^, 
but without the ability of cool silicate dust at 

2 < 10 to soften 7 around ~ 10® cm~^. Evidence 
is mounting that many carbon-enriched metal de- 
ficient stars are in binaries, where they may have 
picked up enriched ejecta of an intermediate-mass 
companion (c.f. Lucatello et al., 2004). The ob- 
served stellar abundance pattern is then still in- 
dicative of the elemental interstellar gas composi- 
tion at early times. 

4.2. Primordial Limit 

In the limit of very low metallicity, Z/Zq < 

3 X 10"'' or less, the work of Schneider et al. (2002, 
2003) and Bromm et al. (2001) suggests that a 
metallicity of ^ 1Q~'^Zq at a density of ~ 10""^^ 
cm~^ leads to efficient cooling and causes a soften- 
ing of the EOS that gives rise to fragmentation and 
the formation of clumps with masses < IMq. We 



would like to stress that the polytropic index 7 is 
already less than unity at much smaller, nu < 10® 
cm~^, densities (Scalo & Biswas 2002) as shown 
by models 1 and 13. While in the density regime 
around 10^^ cm~^, individual clumps are no longer 
able to separate and accretion/coagulation domi- 
nates, at low densities the higher fragment mass 
(Mj oc /0~°'^) may still lead to low-mass fragmen- 
tation (e.g., through disk formation) because the 
medium is very unstable to (turbulent) compres- 
sion when 7 ft: 0.7 (Li et al. 2003). In this case, 
a metallicity of around 3 x 10~'^Zq would still 
retain a soft EOS^, while it would also strongly 
lower the ambient cloud temperature (Bromm & 
Loeb 2003) and reduce the thermal Jeans mass 
(Mj oc T^-^) that pertains to a cloud prior to 
fragmentation. Three-dimensional hydrodynam- 
ical simulations indicate that it is the initial Jeans 
mass that controls, in part, the subsequent frag- 
mentation (Vazquez-Semadeni et al., 1996). Fol- 
lowing the discussion in Spaans &: Silk (2000) on 
the role of the Jeans mass, we would argue that 
7, already at low densities, determines the shape 
of the IMF, while the instantaneous value of the 
Jeans mass during fragmentation plays a lesser 
role. 

4.3. Starbursts 

A similar scenario may apply to starburst 
galaxies. The super-solar abundances (Barthel 
2004), high gas temperatures and high densities 
in those regions imply large amounts of H2O. Also, 
the dust (set here at 75 K) is warm due to the in- 
tense ultraviolet radiation field, which facilitates 
a) H2O heating around 10"* cm~'^ and b) gas-dust 
heating around 10^ cm"*^. Note that H2O col- 
lisional de-excitation dominates the heating for 
super-solar oxygen abundances. In fact, an EOS 
as stiff as in Spaans & Silk (2000) for dust of 100 
K is obtained in this work. 

The resulting IMF, as also argued in Spaans 
&: Silk (2000) for solar abundances, would be 
top-heavy and the overabundance of UV-luminous 
massive stars would continue the starburst cycle 
through radiative feedback. It is interesting to 
note that the starburst model (14) attains a value 

^In a medium with carbon and oxygen in atomic form, fine- 
structure cooling would, in fact, lower 7 to below 0.6 at 
densities below 10* cm~^. 
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of 7 ~ 1-4 at riH ~ 10'* cm~^. This is somewhat 
larger than the canonical value of 7c = 4/3 for 
spherical {v = 3) clouds. 

5. Future Work 

There are a number of issues that should be 
addressed in the future: I) Observations of water 
lines with the HIFI instrument on Hcrschel should 
be performed to determine the H2O level popula- 
tions and substantiate our claim of a polytropic in- 
dex 7 > 1 in starburst galaxies. II) Numerical sim- 
ulations with density and metallicity-dependent 
values of 7, i.e., a pointwise polytropic EOS, 
should be performed for collapsing gas clouds to 
assess the relative impacts of magnetic fields, tur- 
bulence and the EOS on gravitational collapse and 
the shape of the IMF (c.f. Klessen et al. 2004; Li 
et al. 2003). Ill) The difference in influence of an 
external FUV (O and B stars, popIII objects) or 
X-ray (QSOs, mini-quasars) radiation field on the 
EOS of molecular clouds should be studied (Mei- 
jerink & Spaans 2004). 

We would like thank .lan-Willcm Pel for dis- 
cussions on the nucleosynthesis of massive stars 
and Dieter Poelman for help with the excitation 
of water. JS gratefully acknowledges the hospital- 
ity of the Kapteyn Institute where this work was 
performed. 
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Fig. 4. — Figure 4 shows models 1, 14, A, Bl 
and B2; note the progression in 7 as the oxygen 
abundance rises and the dust temperature is larger 
or smaller than the gas temperature. 
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